Superconducting states in the tetrahedral compound PrOs4Sbi2 
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We find possible superconducting states for tetraliedral (T^) symmetry crystals with strong spin- 
orbit coupling using Landau theory. Additional symmetry breaking within the superconducting state 
is considered. We discuss nodes of the gap functions for the different states, secondary superconduct- 
ing order parameters and coupling to the elastic strain. By comparing our results to experiments, we 
find that superconductivity in PrOs4Sbi2 is best described by the three-dimensional representations 
of point group T^. 
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I. INTRODUCTION 

The discovery of superconductivity in the heavy 

fermion compound PrOs4Sbi2^ has spawned a flurry of 
experimental3i4i5i6/,8,9,io,ii,i2,i3,i4,i5,i6,i7,i8,i9 ^^eo- 

retical^i^°i^^i^'^i^3i^^i^''^ activity. PrOs4Sbi2 is the first Pr- 
based heavy fermion superconductor and the first among 
the family of rare-earth filled skutterudite compounds. 
The onset of superconductivity occurs at Td — 1.85 K; 
an additional phase transition is observed as an anomaly 
in the specific heat- and magnetisationiSi at Tc2 = 1.75 K. 
Thermal conductivity measurements in a rotating mag- 
netic field revealed the presence of nodes and a lowering 
of the symmetry of the gap function from four-fold (A- 
phase) to two-fold (B-phase) at T'c2-* Even more intrigu- 
ing is the observation of broken time-reversal symmetry 
in the superconducting (SC) state. By all indications, 
PrOs4Sbi2 is a new kind of unconventional superconduc- 
tor. 

A central issue in the study of unconventional super- 
conductivity is the symmetry of the SC order parameter. 
The phenomenological Landau theory approach is partic- 
ularly useful when little is known about the mechanism 
of superconductivity at a microscopic level, and is ideal 
for describing multiple phase transitions as is the case 
for PrOs4Sbi2. The starting point is knowledge of the 
crystal symmetry group, according to whose representa- 
tions order parameters are classified. The outcome of 
this approach is detailed knowledge of all possible phase 
diagrams and symmetry properties of the SC state, in- 
cluding nodes of the gap function. Phenomenological 
theory can also predict the order of the phase transition. 
While the normal-to-SC phase transition is expected to 
be second order since third-order terms in the Landau po- 
tential expansion2I are prohibited because of gauge sym- 
metry, this is not generally the case for phase transitions 
within the SC state. 

Several theoretical models of the SC order parameter 
in PrOs4Sbi2 have been proposed in order to account for 
the experimental data. Goryo suggested different com- 
binations of s- and d-wave gap functions for the A- and 
B-phases2i in order to account for the change in symme- 
try observed in the thermal conductivity experiment 



The A-phase was assumed to have an anisotropic s-wave 
gap function that has six minima along the [100], [010], 
and [001] directions. In the B-phase, an s + id z2-x2 -way e 
combination was proposed. Different s + g-wave basis 
functions were proposed by Maki et al. for both states.— 
An /-wave pairing state with weak spin-orbit coupling 
was proposed by Ichioka et al. to describe a state with 
point nodes on all three axesJl Finally, Miyake et al. 
considered a microscopic model based on quadrupolar 
fluctuations and nesting in the Fermi surface, and argued 
in favour of Px + ipy-wave pairings^' 

While the models mentioned above may describe par- 
ticular experiments, they can only be considered as em- 
pirical. There are at least two fundamental shortcomings: 
(i) the models are in fact based on the assumption that 
the point group crystal symmetry is Oh- PrOs4Sbi2 has 
lower Tfi symmetry (space group /m3, T^); (ii) there is no 
physical reason why the system should choose one partic- 
ular combination of the basis functions of the irreducible 
representation of the symmetry group over the others. 
Strictly speaking, the theory allows all basis functions 
to contribute to the gap function. Moreover, the coeffi- 
cients in the such combinations in general depend on the 
external conditions (temperature, magnetic field, etc.). 
Only such a general state is thcrmodynamically stable 
and occupies a finite region of the phase diagram. 

In this paper we use the Landau theory approach to 
classify SC phases for tetrahedral (Th) crystals, includ- 
ing those which may be reached by additional symmetry 
breaking within the SC state. We use the strong spin- 
orbit coupling limit in which the spin rotation symmetry 
is broken. 

The first attempt to accomplish such a classification 
was made by Gufan;^ In Sec. II of this paper we use a 
different approach and reproduce most results of Ref. 113 
for Th symmetry— In addition, we discuss the basis func- 
tions of the irreducible representations, the gap function 
nodes and the orders of the phase transitions between 
different SC states. In Sec. Ill we consider secondary SC 
order parameters which influence the nodes of the gap 
functions. In Sec. IV the coupling between the SC order 
parameters and elastic strain is discussed. Sec. V is de- 
voted to matching the experimental data with the states 
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found theoretically. Sec. VI summarizes the paper. 



II. CLASSIFICATION OF 
SUPERCONDUCTING STATES 

A procedure for constructing SC classes and finding the 
gap nodes with strong spin-orbit coupling was originally 
proposed by Volovik and Gor'kov (VG)r& who listed all 
SC states which can be reached from the normal state 
by a second order phase transition for Oh, -D4/1 and Dq^ 
crystals. One begins by classifying possible order pa- 
rameters according to the representations of the crystal 
point group. In systems with inversion symmetry, all 
representations have a definite parity. Those with even 
parity must be matched with singlet pairing of the spin 
states for the pair wavefunction to be antisymmetric; like- 
wise odd parity representations are matched with triplet 
spin states. For each parity, the group Th has a one 
dimensional representation A, a two dimensional repre- 
sentation E, which is reducible to two one-dimensional 
representations that are complex conjugate, and a three 
dimensional representation Tm^ 

The SC gap function is a 2 x 2 matrix in pseu- 
dospin space given by A(k) = iayipih) for singlet pair- 
ing and by A(k) — i{d{'k)a')ay for triplet pairing, where 
S' = ($x,^y,^z) are Pauli matrices, V'(k) is an even 
scalar function and ci(k) is an odd pseudovector func- 
tion, ^(k) and d(k) are expressed in terms of the com- 
ponents of the order parameter rji as ipi}^) — J2i Vi^iO^) 
and d(k) = '^■r]idi{'k). Here ipi(k) and di(k) are the 
basis functions for the even (spin-singlet case) and odd 
(spin-triplet case) irreducible representations of the point 
group, respectively 

The method of finding the SC states implemented by 
VG is to construct a Landau energy functional of rji for 
each order parameter that is invariant under G x U x /C, 
where G is the point group, U is gauge symmetry and 
K, is time-reversal, and analyse its extrema. In order 
to account for all possible phase diagrams, a very large 
number of terms must be included and the analysis of 
such a cumbersome model is tedious at best.'^^ In prac- 
tice, terms are restricted to those needed to describe the 
normal to superconducting phase transition(2Si^ while 
states resulting from additional phase transitions within 
the SC state must be found by other methods. The VG 
approach can be applied Th crystals. However, here we 
use even an easier method. We use the fact that Th is a 
subgroup of Oh- Therefore we start with the results for 
the symmetry groups of SC classes obtained by VG for 
Oh symmetry and reduce them by removing the symme- 
try elements that are absent in the normal state of Th 
symmetry. 

We consider additional symmetry breaking within the 
SC state by constructing effective Landau functionals 
of effective order parameters, which describe the phase 
transitions between SC states with a group-subgroup re- 
lation. This procedure is straightforward, since the sym- 



metry group of a SC state is discrete (the continuous 
gauge symmetry is already broken). In the following, 
we consider the two dimensional representation in detail, 
while only the results are given for the three dimensional 
representation. 

Our results are summarized in Table ^ It lists all 
possible SC states for both even parity and odd parity 
which are possible if only a single irreducible represen- 
tation is present. We define the relations between the 
components of the order parameters, the symmetry of 
the SC state, and the structure of nodes in the gap func- 
tion. We make the distinction between accidental, ap- 
proximate and rigorous nodes. Accidental nodes occur in 
empirical models when a particular form of the gap func- 
tion is chosen a priori, such as that proposed in Ref. 12^ . 
Such nodes cannot be stable because even small contribu- 
tions of functions with the same symmetry remove them 
immediately;^^ Accidental nodes are unphysical and so 
we disregard them. Approximate nodes are a property 
of all possible basis functions which can be constructed 
for a given representation. These nodes may be removed 
when admixtures of other representations, which couple 
to the SC state as secondary order parameters, are taken 
into account, thus leaving only rigorous nodes required 
by the symmetry of the SC stateiSS*^ The secondary or- 
der parameters are proportional to the third power of 
the primary order parameter]^ Hence, the experiments 
that probe the symmetry of the gap function close to Tc 
may find the approximate nodes, while only the rigorous 
nodes remain when T — > 0. A more detailed discussion 
of the secondary order parameters is given in Sec. III. 



A. ID representation Ag^u 

The analysis of the one-dimensional representations Ag 
and Au is straightforward. Only gauge symmetry is bro- 
ken and there are no nodes. The symmetry of the SC 
state is T X AC. In the lowest order in k, the basis function 
for the singlet channel V'(k) is constant on the Fermi sur- 
face and for the triplet channel d(k) fc^x + kyy + fc^z. 
Here and below means 'transforms like' so that all 
our results remain valid for higher order basis functions. 



B. 2D representation Eg^u 

We choose the basis functions of the two-dimensional 
representations Eg and E^ in complex form as in Ref. |2^ 

-01 ~ fc^ + eky + e^fc^, ip2 ~ kl -\- e'^ky + efcf ; 
di ~ fcxX + skyY + e^fc^z, d2 ^ fc^x -f s^kyy -f efc^z, 

(1) 

where e = exp(27ri/3). With this choice of the basis 
functions, the order parameter has the following trans- 
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TABLE I: SC states described by one irreducible representation of the point group T^. The relative magnitudes and phases 
of the components of the order parameter are defined in the first column. The symmetry groups of the SC states are listed 
in the second column. Approximate and rigorous nodes of the gap function for even parity are listed in the third and fourth 
columns, similarly for odd parity in the fifth and sixth columns. The word 'same' is used when rigorous nodes coincide with 



approximate nodes. The numbers 




parentheses in the fifth column indicate whether 


onG or 


UVJLll tiCLUO IICLVC lHJvJ-Co. 




State 


Symmetry 




Approximate 


Rigorous 




Approximate 


Rigorous 








nodes 


nodes 




nodes 


nodes 


(1) 


TxlC 




none 


none 




none 


none 


(1, 0) 


T{D2) 




8 points (111) 


same 




8 points (111) (1) 


same 


(</>!, <f>2) 


DiXJC 


E, 


8 points (Iff) 


none 


Eu 


none 


none 




D2 




8 points (111) 


none 




none 


none 


(1, 0, 0) 


D2{C2) X /C 




2 lines ky = 0, = 


same 




2 points [100] (2) 


same 


(1, 1, 1) 


Ca X /C 




6 points (001) 


none 




none 


none 


(1, e, e^) 


C3{E) 




6 points (001), 2 points [111] 


2 points [111] 




2 points [111](1) 


same 




D2{E) 




1 line fc^ = 0, 2 points [001] 


same 




none 


none 


i\ni\,\r)2m 


C2{E) X /C 




1 line k^ =0,2 points [001] 


same 


Tu 


none 


none 




C2{E) 




1 line k^ =0,2 points [001] 


same 




none 


none 


(l»7iMI'?2|, |r73|) 


C'2{E) 




6 points (001) 


none 




none 


none 


(l^il> M, \m\) 


K. 




6 points (001) 


none 




none 


none 




E 




6 points (001) 


none 




none 


none 



formation properties, 

U{0){rii,V2) = e^\m,m). (2) 

where C2 stands for any of the twofold rotations in Th, 
and Cg^^ is a 27r/3 rotation about the [111] direction, 
and U (9) is a gauge transformation. 

In Table I, three states are listed for the two- 
dimensional representations of Th ■ These differ from the 
Oh states (1,0), (1, 1) and (1, —1). As shown below, the 
extra freedom in the phase and magnitude of the last two 
states of Tfi arise from terms in the free energy which are 
allowed under Th but not Oh- 

The SC state (1,0) in Oh corresponds to the group — 

0{D2) = {I?2,2C|/C,2C|i7(27r/3)/C,2C|t/(4V3)/C, 
2Cf /C, 2C72^^C/(27r/3)K:, 2C|^[/(47r/3)/C, 
4C3C/(M3),4Cff/(2V3)}, (3) 

where D2 is the group of two-fold rotations about the 
[100], [010], and [001] axes. In Th, the remaining symme- 
try elements are 

T{D2) = {D2, 4[/(47r/3)C3, 4f/(27r/3)C2}. (4) 

Considering the symmetry groups of the states (1,1) 
and (1, —1) in Oh, which are D4 x /C and 04(02) x K, 
respectivelyfS& where 

04(02) X /C = {O2, 2C^U(n), 2Cf [/(tt)} x JC, (5) 



we notice that they both reduce to the same symmetry 
O2 X K, in Th- Moreover, it follows from ||2J) that this 
symmetry does not fix the relation between the phases 
and (1)2 of the OP components 771.2 = 1^71,21 exp(i(/)i_2), 
but the magnitudes are equal — \ri2\. Therefore, 
we denote this state as (4>i,4>2)- This also may be ver- 
ified from the following Landau model which describes 
the Eg_u representation of Th'. 

F = a(\m\' + \V2\') + l3i(\vi\^ + \v2\') + 2f32\in\'\v2\' 

+ii{vhf + vhf) + i2i(vfvf - vhf), (6) 

where a, /32, 71, and 72 are phenomenological parame- 
ters. The last two terms reduce to 271 [771 ['^ [772 1^ cos(3(^) -I- 
272|'7iP|?72|'^ sin(3(/)), where (j) = — (j)2- Thus the equi- 
librium value of (j) depends on the (generally temperature 
dependent) ratio 71/72. In contrast, in Oh symmetry the 
72 term is prohibited, hence (f) is fixed to be either [for 
(1, 1) state] or tt [for (1, —1) state]. 

The gap function of the (0i,02) state in the singlet 
channel is^i 

A(k) - cos((/)/2)A:2 +cos(0/2 + 27r/3)A:^ 

+ cos(<j>/2 + 4TT/3)kl (7) 

and in the triplet channel, 

A(k) - cos^((l3/2)kl+cos^((t)/2 + 2TT/3)kl 

+ cos^(4>/2 + ATT/3)kl. (8) 

We would like to stress that the state dx2-y2 and its 
equivalents, obtained by permutations of x, y, and z, are 
not stable in Th- Instead, they are replaced by the more 
general state (0i, 02 ) with the gap function {T)). 
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FIG. 1: Second order phase transitions among states of the 
Eg and i5„ representations of T^. 

Normal state 

(1,0) (01,02) 

In Oh, the states (1,0), (1,1), and (1,-1) are 
connected to the normal state by a second order 
phase transition^ Since up to fourth order terms, the 
model (jnji coincides with that of Oh, we conclude that 
the states (1,0) and (</>i,</>2) can be reached from the 
normal state in Th by a second order phase transition. 

There is a third state which can be described by the 
E representation in Th- Its symmetry group is D2 (time 
reversal is broken) , which is a common subgroup of both 
T{D2) and D2 x IC. As is seen from it has no con- 
straints on either the magnitudes or phases, therefore we 
denote this state (771,772). 

The phase transitions (1,0) (771,772) and (0i,(/)2) 
(7^1,772) may be described by effective free energy func- 
tionals of effective order parameters. The factor group 
T{D2)I D2 is isomorphic to the cyclic group C3 and thus 
third order terms are allowed in the effective free energy 
describing the (1,0) (771,772) transition. This implies 
that it cannot be second order. The effective free energy 
is a functional of the effective order parameter 7^2, 

Fe^/[(1,0)^ (771,772)] = 5|772p + 71(77! +77*3) 

+^72(ry|-r/2*')+/3|772|'(9) 

On the other hand, a second order transition (0i, ^2) — > 
(?7ij?72) is possible. This transition is described by an 
effective order parameter 6 = [771 1 — I772I. Third order 
terms are prohibited in the effective free energy because 
of time reversal symmetry, 

F,ffU^,cj,2) ^ (7/1,772)] = a'5^^p'5\ (10) 

There are no other states described by the E represen- 
tation alone, because the basis functions ^ are invariant 
with respect to all symmetry operations of D2 group and 
there are no other symmetry groups containing D2. 



C. 3D representation Tg^u 

The lowest order basis functions for the Tg representa- 
tion of Th are 'd-wave' (i.e. second order in k), 

■01 ^ kyk^, -02 k^k^, 03 kxky, (11) 

while for the representation the lowest order basis 
functions are 'p-wave', and there are two independent 



sets of them, 

di ~ akyZ + bk^y, 

d2 ^ akzSc + bkxZ, 

ds ~ ak^y + bkyX. (12) 

Here a and b are arbitrary numbers, in contrast to Oh, 
which fixes b — —a in Tiu representation and b — a in 
T2U representation. It follows that the order parameter 
transforms as 

C2{Vl,V2,V3) = {-Vi,-V2,V3), 

C3"(77i, 772,773) = (772,773,7/1), 

/C(77i, 772,773) = (771, 7/*, 7/*), 

U{9){m,V2,V3) = 6^^(771,772,7/3). (13) 

To find the SC states of the three dimensional repre- 
sentation, we again use the Oh states as a starting point. 
For Oh, there are four states accessible by a second order 
phase transition from the normal state: (1, 0, 0), (1, i, 0), 
(1; 1, 1) and (1, £, s'^), with symmetries 

D^iCi) X /C = {E,C^,2C^,AU{tt)C^''} x IC, 
Di{E) = {E,U{t:)C^,2U{±tv/2)CI,C^1C, 
C/(^)C|/C, 2U{±TT 12)02" JC}, 
DsiCs) X /C = {E,2C3,3U{tt)C2^} x /C 

DsiE) = {E,U{A7:/3)C3,U{27r/3)Ci,CrJC, 

C/(27r/3)C2"^^/C,C/(47r/3)C2^/C}, (14) 

respectively Here E is the identity element. Reducing 
these groups, we find the following classes for Th'- 

02(02) xJC = {E,C'^,U [1^)01, U{i:)C^} X K. 

D2{E) = {E,U{7r)Ci,C^IC,U{7r)ClJC} 

C3 X /C = {E,C3,Ci} X IC 

CsiE) = {E,U{A7r/3)C3,U{27r/3)Cl} (15) 

We notice that the D2{E) symmetry actually does not 
require |77i| = I772I. Hence, the state (l,z,0) is not stable 
in Th- Instead, it is replaced by the state (I7711, 7|772|, 0). 
A direct second order normal-to-(|77i|, i|7/2|, 0) transition 
is possible in Th- These findings are also evident in the 
form of the Landau potential for the 3D order parameter. 
In order to display the Th (but not Oh) symmetry, a 
Landau model for Tg and T„ must include at least sixth 
order terms, as in the case of Eg and Eu- These sixth 
order terms are composed of five linearly independent 
invariants: 

\vi\' + \mf + \v3\', \m?\v2\'\m\', 

(Imp + Im? + \v3\''){v1vf + vh3'' + vlvf + c.c.), 

i\Vi\'\V2? + \ri2\'\ri3\' + \V3\'\vi\') 

±{\mW + \v2?\v3\' + \v3\'\m\'), 

ivtvf + Vhf + Vhf) ± (77|77f + vhf + vhf) + c.c. 

(16) 

The negative signs in the last two invariants in Hlt)|) occur 
in Th but not in Oh- 



5 



FIG. 2: Second order phase transitions among states of the 
Tg and r„ representations of T^. 

Normal state 
(1, 0, 0) (1, 1, i)\^(i> f') 



i\n\\, \ri2\, Iml) Ui, 0) /|r72|, l'?3l) 

(f7b Ui^ 13) 



Considering all possible subgroups of the groups 
in H15|l . we find five more SC states as listed in Tabled 
where 



C2iE) 

cm 



{E,Ui{^)CllC}. 



(17) 



We have examined the transitions within the SC state by 
considering effective free energies which describe them, 
similar to those described for the 2D order parameter, 
Eqs. lO, IjlOl) . The diagram of all second order phase 
transitions described by the three-dimensional represen- 
tations of Th is given in Fig. |21 

Thus we find that the absence of four-fold rotation 
symmetry in PrOs4Sbi2 essentially changes the struc- 
ture of possible SC states. The states (1, 1) and (1, —1) 
are not stable, because the value of in Eq. {7)) is not 
fixed. Similarly, the state (l,i,0) is absent in the three- 
dimensional representations. Additionally, all SC states 
which may be connected to the normal state in Oh, D^h 
or Z?g/i symmetry by a second order phase transition are 
one-parameter in the sense that all components of the 
order parameter are proportional to one quantity, its ab- 
solute valuei2SiM The situation is different for the states 
{4>ii4>2) and (jryil, i|772|, 0) in T/j, for which two indepen- 
dent quantities describe the SC state. 



III. SECONDARY SC ORDER PARAMETERS 

In general, the primary order parameter is accompa- 
nied by secondary order parameters which do not change 
the symmetry of the SC state. The influence of secondary 
order parameters on the gap nodes was discussed in Sec. 
I. Since secondary order parameters do not change the 
overall symmetry of the superconducting state, they are 
most easily found by identifying supergroups of the states 
listed in the second column of Table ^ which correspond 
to another superconducting state. Table ITTl lists them. 

In order to calculate how the secondary order param- 
eters appear in the ordered phases, we need invariants 
of the types rp^ and rf^"^, where ry is the primary or- 
der parameter and ^ is the secondary order parameter. 
From the first type of invariant, it is clear that ^ and r] 



must have the same parity. There are three scenarios to 
consider: i) The 2D primary order parameter with ID 
secondary OP ii) 3D primary with ID secondary and iii) 
3D primary with 2D secondary. In the rest of this section, 
we denote the primary order parameters as rjj — \r]j\e'^'^^ 
and the secondary order parameters as = 10 1^'^^- 



A. 2D primary with ID secondary 

The coupling terms of the two order parameters in the 
Landau potential are 



-irifv2C + V*ivU*)], 

{\vi\' + \v2m?- 



(18) 



In the state (1,0), the first two terms vanish, hence ^ = 
0. In the state ((/>i,(^2), the first two terms are finite 
and 1^1 cx \r]\'^ . Minimization with respect to 6 yields 
9 = 5(01 + 02)- This relation between the phases of the 
OP's ensures that time reversal symmetry in preserved. 
There is no such relation between the phases when the 
primary order parameter state is (?7i,?72). This reflects 
the fact that time reversal symmetry is broken. 



B. 3D primary and ID secondary 



The coupling terms are 



{mv2V3 + vlmvs + vlnlvsi^. + cc, 

{vf + r]f + rif)e + cc. 



(19) 



(hil 



It follows that if any of the components of the 3D order 
parameter is zero then the potential has a minimum at 
^ = 0. This is also the case for the state (l,e,e^). In 
the states in which (pi = (j)2 = 4>s [*• e- (IjIjI) and 
(|'7iM'?2|, |?73|)], one obtains = 0i. However, in the 
state (|?7iMI'72|, byal) we find 6* = 0i ± f . 



C. 3D primary and 2D secondary 



The coupling terms are 

(?7i*'72'73 + e?7i'72^3 + £^??i?72?73)Ci* 

+iviV2V3 + svlm-nl + £^?h* ^2^3)6 
i[(jilmm + £Viii2V3 + e'^mrnvDCi 

-{V1V2V3 + £v*iV2V3 + s'^Ti*i-nlm)i2 
ivf + vf + %*')ei*a + cc, 

{\m?+e\v2?+s'\V3?)iiC2+c.c, 

\m? + \v3\'mi? + \^2\ 



cc. 



(20) 



{\m? + \m? + \v3\'mi? + \^2\'). 

For this type of mixing we only consider the (l,e,£^) 
state of the primary order parameter, since in the other 
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TABLE II: Secondary SC order parameters. The primary SC 
order parameters are listed in the first column and all sec- 
ondary SC order parameters are listed in the second column. 



Primary 


Secondary 


(1) 




none 


(1,0) 




none 


(01, 02) 




(1) 






(1) 


(1,0,0) 




none 


(1, 1, 1) 

V -^7 -^7 




(1) 


(!,£,£') 




(1,0) 


(l'7iMI'72|,0) 




none 


i\vi\,\m\,o) 




none 






none 


ilml-Mvil, \v3\) 


(1) 


, (01,02) 


(Iml, M, M) 


(1) 


, (01,02) 




(1) 


, ('7l,'72) 



A. ID order parameter 

There is no difference between Oh and in this case. 
The coupling of the strain to the SC order parameter is 
described by the following term in the Landau potential 

Frje^ p\v\^{ei + e2 + e3) (22) 

The dilatation al strain 61+62+63 appears as a secondary 
order parameter, and the only elastic constant which is 
discontinuous is Cn. 



B. 2D order parameter 

The coupling terms are 

Frje - Pl(|»7lP + l'72n(ei +62 +63) 

+P2[»7lf?2(ei + £62 + £^63) + C.C.] 

+ip3[viV2{'^i + £62 + £^63) - e.c.]. (23) 



states where Eg^u is present as a secondary order param- 
eter, „ is also present, and it surely removes all nodes. 
The first two invariants in the (1, e, e^) state reduce to 
6I771PIC2I cos(6'2 -0i) and 6|?7i|^|^2| sin(6'2 -0i), respec- 
tively. Thus, the state (0, 1), which is equivalent to (1, 0), 
appears as a secondary effect. Note that 62 — (f>i is not 
fixed, which is expected since the state breaks time re- 
versal symmetry. 



The third term is absent in 0^- The free energy of the 
OP is given by Eq. 10, which describes the second or- 
der phase transitions between the normal state and the 
superconducting states (1,0) and (0i,(/)2)- 

Diatomic strains 62 — 63 and 2ei — 62 — 63 appear in 
the transition to (01,02)- Therefore, it is necessary to 
average the elastic moduli in all three directions to take 
into account domains. 



IV. STRAINS AND ELASTIC MODULI 

Unconventional SC states normally break spatial sym- 
metry in addition to gauge. If the crystallographic class 
changes, one can expect the development of new com- 
ponents of the strain tensor and certain anomalies in 
the elastic moduli which can be measured by ultrasound 
propagation.^* Such a measurement has not yet been re- 
ported for PrOs4Sbi2. Thus, here we consider all repre- 
sentations for the normal-to-A phase transition. 

The elastic energy for Th is the same as for Oh, 

Fel - ^(e?+6^+62)+Ci°2(eie2+e263 + 6i63) 

+ %el + el + el), (21) 



C. 3D order parameter 

The coupling terms are 

F,^e = Pl(hl|^+l?72|^ + |'73|^)(ei+62+e3) 

+P2[3(h2p-h3n(e2-e3) 

+ (2|ryi |2 - |,y2 I' -|?73 I') (261 - 62 - 63)] 
+P3[(|??2|^ - |?73|^)(26i - 62 - 63) 

-(2|,?i|2-|,/2|'-|%n(e2-63)] 

+P4[{V2V3 + mvl)ei + ijilvi + v^.-nDe-b 
Hvlm + mv2)E&]- (24) 

The third term appears in Th but not Oh- Shear strains 
64,5,6, but not diatomic strains, are present when all three 
components of the OP have the same magnitude. Di- 
atomic strains appear when any of the magnitudes differ. 



where 61,. ..,6 are the components of the strain. Generally, 
if the strain is a secondary order parameter it couples 
to the primary order parameter as rf'e, which leads to a 
development of the secondary order parameter as 6 ~ 77^ . 

The development of the strains following each normal- 
to-SC transition and discontinuities of the elastic moduli 
are shown in Tabic UTTI 



V. DISCUSSION 

Experimentally, the symmetry of the SC states and 
the nature of the phase transition between them in 
PrOs4Sbi2 is far from resolved. Anomalies at Tc2 have 
been observed in many experimenta^iSiLSiiSiiiiSiiSiiS. 

Specific heat measurements by Volume et al& found a 
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TABLE III: Strains and discontinuities in the elastic moduli following normal-to-SC phase transitions in Th crystals. The SC 
states are listed in the first column. Strains which appear as secondary order parameters and discontinuities of the elastic moduli 
are listed in the second and third columns respectively, as functions of the primary order parameter and the phenomenological 
constants. The fourth-order coefficients Pi in the Landau potential for the 2D order parameter a defined in Eq. JHJ. For the 
ID and 3D order parameter they correspond to the following terms^ /SI*?!* and /3i(|7;i|^ + \ri2'[^ + |?73|^)'^ + /32|77i + ?72 + Jyil^ + 
/33(|r;ip|T?2p + |?72p|f?3p + l^yiPl'/sP); respectively. The domain average values for the elastic moduli dj are calculated as 
C^V = {Cu + C22 + C33)/3, C^2 = {C12 + C23 + Ci3)/3. The superscript denotes the values in the normal state. 



Transition: 


Strains which appear 


Elastic moduli 


Normal to 


as secondary order parameters 


in the SC state 


(1) 


1 1 — 3p T) 1 

ei + 62 + ea = 0^^+20"^ 


Cll = C22 = C33 = C*!! — ^ 

Cii — C12, C44 continuous 


(1,0) 


ei + 62 + 63 = r'^''2l% 


Cll = C22 ~ C33 = Cll — 
Cll — C12, C44 continuous 



r» —6|r?i 1^(09 cos (h — pr, sin d>') 

261 - 62 - 63 = — ' ^wi — -Zr^ ^ 

sin 0+P3 cos (j)) 



62 — 63 



11 — L^ll 



c 



Av 



12 - ^12 - 2{20i+fe) 

C44 continuous 



'-■12 



(1,0,0) 



61 + 62 + 63 = 



261 

62 - 



62-63 



-3pi|i7i 
'll+^'-'12 
-12p2l'7ll^ 



CV 



^3 r-c^Lr"'" 



3pj + 24p^+8p^ 
6(/3i+/32) 
_ 3pj-12p|-4pi 
12 ^12 6(/3i+/32) 

C44 continuous 



>-^ll 

c 



L-ii 

^12 



(1,1,1) 



61 + 62 + 63 = ^1^^ 



64,5,6 



2P4l'7ir 
' — r" — 



Cll = C22 ~ C33 = Cll 
Cll — C12 continuous 

1^ _ f~iO 2p4 

044 - O44 3(3/3i+3fe+/33) 



JpT 



2(3/3i+3/32+/33) 



(l,£,e') 



61 + 62 + 63 = 



-9pi\vi\' 



Cll = C22 ~ C33 = Cll 

Cll — C12 continuous 



2(3/3i+/33) 



C44 = C4, 



P4 



6(3/3i+|33) 



(l'7lMI'72|,0) 



61 + 62 + 63 = 
61 + 62 - 263 
61 - 62 



-3pi(t7;^+i7S) 
Cjj+2Cj2 

-6P2(I'71P + I'72l^)+6P3(I'?1| 



7^ 



i')2r) 



3pj(4/32-/33)+4(3p^+p^)(6/3i+2/32+|33) 

6(4/32 -ft )(4/3i+/33) 

-^Av _ ^0 3pf(4/32-/33)-2(3pi+pi)(6/3i+2/32+fe) 
^12 — (--12 



6{4/32-/33)(4/3i+/33) 



6P2(l'?ll^-]'72i^)-2p3(l'?ll^ + l'?2r) 
Oil — , o 



C44 continuous 



jump at Tc2, indicative of a second-order phase transi- 
tion. On the other hand, Aoki et almii^ found a kink, re- 
sulting in a steeper temperature dependence below Tc2, 
which seems to correspond to a first-order phase tran- 
sition. The most dramatic observation is the change in 
symmetry at the A-B phase transition seen in thermal 
conductivity measurements.^ The double transition was 
also observed in magnetisation measurements as a peak 
effect in M{H)m^ii^ One of these measurements found 
strong anisotropiesfi^ possibly indicative of a change in 
symmetry; the other did not^ Finally, recently pene- 
tration depth measurements have been interpreted not 
as a phase transition, but rather as a cross-over due to 
two-band superconductivity,^^ 

The temperature range in which the A-phase exists 
is very narrow, thus with two exceptions jSii^ the re- 
ported experiments probe the properties of the gap in 
the B-phase. Experiments consistently rule out the ex- 
istence of line nodes in the B-phase However, 
the presence of point nodes in the B-phase is clearly in- 



dicated by a power law temperature dependence of the 
specific heatj2i the thermal conductivity measurement)^ 
and the penetration depthfiS. Nuclear quadrupolar res- 
onance experiments'^ can be interpreted as either fully 
gapped or nodes. Very low temperature tunneling spec- 
troscopy-^^ finds no nodes at all in the B-phase, perhaps 
consistent with rigorous nodes, rather than approximate 
nodes. Finally, /iSR^ indicates that the B-phase is fully 
gapped. 

Only a couple of experiments have specifically dealt 
with the symmetry of the gap function. ^'^^ In the ther- 
mal conductivity experiment, point nodes were found in 
the [010] direction in the B-phase and in both the [100] 
and [010] directions in the A-phase.^ However, in this 
measurement, there is no clear explanation for why the 
two-fold symmetry is actually observed as such, rather 
than averaged out into domains. The penetration depth 
has a power law temperature dependence corresponding 
to point nodes along all three principal crystallographic 
axes»i£ No studies of the nodal structure along the [111] 
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direction have been reported so far. An extremely impor- 
tant finding is due to another /xSR measurement, which 
showed that time reversal symmetry is broken in the B- 
phaseii^ 

In determining which of the states listed in Table I best 
describes PrOs4Sbi2 we make the following assumptions: 

i) the B-phase breaks time reversal symmetry; ii) there 
are point nodes in the B-phase located in the [100] and/or 
equivalent directions, and there are no line nodes in the 
B-phase; in) the A-B phase transition is second order; iv) 
both phases are described by the same order parameter. 
The first two assumptions are based on fairly conservative 
interpretations of the experimental data available to date. 
We use the last two assumptions to narrow the choices 
of possible states. Their validity is subject to further 
experimental study. 

We exclude the A and E representations because of 

ii) . In Tg and T„ representations, the first four states 
listed in Table are connected to the normal state by 
a second-order phase transition (see Fig. I^J , but among 
them only (1,0,0) and (Ir^il, z|?72|, 0) may be followed by 
another second order phase transition involving the same 
order parameter. Therefore, these are the only two pos- 
sibilities for the A-phase. If the A-phase is (|?7i|, i|?72|, 0) 
then the B-phase is either (771,772,0) or (|r/i|, i|772|, lr/3). 
The former is excluded because it has line nodes in the 
singlet channel and no nodes at all in the triplet chan- 
nel. The latter possibility must be singlet because it has 
no nodes at all in the triplet channel. If the A-phase is 
(1,0,0) then the B-phase is the (|7;i|, 7|772|, 0) state. Be- 
cause there are no line nodes in the B-phase, the pairing 
is therefore triplet. Strictly speaking, (|77i|,7|7;2|,0) has 
no nodes at all under symmetry. However, nodes ap- 
pear in the corresponding Oh state (1, i, 0)^ Such nodes 
may be pronounced dips in Th if the Fermi surface has 
the approximate Oh symmetry as found in Ref . |3 There- 
fore, the two most likely possibilities for the sequence of 
SC phase transitions in PrOs4Sbi2 are 

normal {\vi\,i\V2\,0) ^ {\Vi\,'i\V2\, 1%!) 

in the singlet channel and 

normal (1,0,0) -> (|r7i |, 7|772|, 0) 



in the triplet channel. 

If future experiments fail to be consistently described 
within the framework described in this article, then it is 
likely that the assumption that the order parameters of 
both transitions belong to the same representation will 
merit closer examination. It is possible that the B-phase 
may be due to the appearance of an order parameter 
that belongs to a different representation than that of 
the A-phase. This possibility is somewhat unsatisfactory 
in situations when the phase transitions occur very close 
together, as in PrOs4Sbi2, because it suggests a rather 
fine tuning of the phenomenological parameters. Second- 
order phase transitions between any states which are re- 
lated as group-subgroups are allowed, provided third or- 
der terms of the effective order parameter are absent in 
the free energy. The order parameter of the B-phase may 
be a superconducting order parameter that belongs to a 
different representation than that of the A-phase, or it 
could even be something completely different, such as a 
structural order parameter, or a state with broken trans- 
lational symmetry. 



VI. SUMMARY 

To summarise, we find group-theoretically the SC 
states which can be realized in crystals with Th sym- 
metry. Additional symmetry breaking within the SC 
state is considered. Heavy fermion superconductivity 
in PrOs4Sbi2 is best described by the three-dimensional 
representations of Th point group. Considering experi- 
mental results, we propose the two most likely scenarios 
for the SC phase transition sequence found in PrOs4Sbi2 , 
one in the singlet and another in the triplet channel. 
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